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Abstract 

The Conte-Musette method has been modified for the search of only elhptic solu- 
tions to systems of differential equations. A key idea of this a priory restriction is to 
simplify calculations by means of the use of a few Laurent series solutions instead 
of one and the use of the residue theorem. The application of our approach to the 
quintic complex one-dimensional Ginzburg-Landau equation (CGLE5) allows us to 
find elliptic solutions in the wave form. We also find restrictions on coefficients, 
which are necessary conditions for the existence of elliptic solutions for the CGLE5. 
Using the investigation of the CGLE5 as an example, we demonstrate that to find 
elliptic solutions the analysis of a system of differential equations is more preferable 
than the analysis of the equivalent single differential equation. 
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1 Introduction 



At present methods for construction of special solutions of nonintegrable sys- 
tems in terms of elementary (more precisely, degenerated elliptic) and elliptic 
functions are actively developed [2,7,8,9,16,17,21,25,26,27,29,31,32,33,37,40,41,43,49] 
(see also [30] and references therein). Some of these methods are intended for 
the search for elliptic solutions only [26], others allow to find either solu- 
tions in terms of elementary functions only [7,8,31,49] or both types of solu- 
tions [2,9,16,17,21,25,27,29,32,33,37,40,41,43]. Note that the methods [40,29] 
allow to find multivalued solutions as well. 
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Elliptic and degenerate elliptic functions arc single- valued functions, therefore, 
such solutions of a nonintegrablc system exist only if there exist the Laurent 
series solutions of it. Such local solutions can be constructed by means of the 
Ablowitz-Ramani-Segur algorithm of the Painleve test [1] (see also [24,30,39]). 
Moreover for a wide class of dynamical systems using this method one can find 
all possible Laurent series expansions of solutions. In this way one obtains 
solutions only as formal series, that is sufficient, because really only a finite 
number of coefficients of these series is used. Examples of construction of such 
solutions are given in [31,42]. The Laurent series solutions give the information 
about the global behavior of differential systems and assist to look for both 
exact solutions [32] and the first integrals [28]. The Laurent series solutions 
can be used to prove the nonexistence of elliptic solutions [23,47] as well. 



In [32] R. Conte and M. Musette have proposed a new method for construction 
of single- valued special solutions of nonintegrablc differential equations. A key 
idea of this method is the use of the Laurent series solutions to transform 
the initial differential equation into a nonlinear system of algebraic equations. 
Using this method one can in principal find all elliptic and degenerate el- 
liptic solutions. Unfortunately if the initial differential equation includes the 
large number of numeric parameters, then it is difficult to solve the obtained 
nonlinear system of algebraic equations. 



The goal of this paper is to propose a modification of the Conte-Musette 
method, which allows to seek elliptic solutions only. We show that in this 
case it is possible to fix some parameters of the initial differential system and 
therefore simplify the resulting system of algebraic equations. To do this we 
use the Hone's method, which has been proposed to prove the non-existence of 
elliptic solutions [23]. Note that using our approach one can find in principal 
all elliptic solutions. 



In [32] the authors have transformed the initial system of two coupled ordi- 
nary differential equations into the equivalent single differential equation and 
only after this have constructed the Laurent series solutions. In our paper we 
demonstrate that the analysis of the system of differential equations may be 
more useful than the analysis of the equivalent differential equation. Moreover 
in this paper we show that if a system of differential equations includes a few 
functions it is possible to find the analytic form of a function which satisfies 
this system even without knowledge of other functions in the analytic form 
and without elimination of them from the system. 
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2 The Conte-Musette method for the system of differential equa- 
tions 

In [32] R. Conte and M. Musette have proposed a way for searching of ellip- 
tic and degenerate elliptic solutions to a polynomial autonomous differential 
equation. In this section we reformulate this method for a system of such 
equations: 

Fi(yl:\ylr'\ y;t, y) = 0, i^l,...,N, (1) 

where y = {yi{t) , y2{t) , . . . , |/L(t)} and y^':} = 

It is known that any elliptic function (including degenerate one) is a solution 
of some first order polynomial autonomous differential equation. The classical 
results of P. Painleve, L. von Fuchs, C.A.A. Briot and J.-C. Bouquet allow 
one to construct the suitable form of such an equation, whose general solution 
is a meromorphic function with poles of order p (see details in [32]) ^ : 

m {p+l){m-k)/p 

Y: E Wyt' = o, ho,m^i, (2) 

fe=0 j=0 

in which m is a positive integer number and hj ^ are constants to be deter- 
mined. The general solution of (2) is either an elliptic function, or a rational 
function of e'^^, 7 being some constant, or a rational function of x. Note that 
the third case is a degeneracy of the second one, which in turn is a degenerate 
case of the first one. 

The Conte-Musette algorithm is the following [32] (see also [9]): 

(1) Choose a positive integer number m, define the form of Eq. (2) and 
calculate the number of unknown coefficients hj^k- 

(2) Construct solutions of system (1) in the form of the Laurent series. If 
such solutions do not exist or they correspond to known exact solutions, 
then no unknown single- valued solutions exist. Note that, since system 
(1) is autonomous, the coefficients of the Laurent series do not depend 
on the position of the singular point. They may depend on values of 
the numerical parameters of (1). In addition, some of these coefficients 
(the number of which is less than the order of system (1)) may take 
arbitrary values and have to be considered as new numerical parameters. 
One should compute more coefficients of the Laurent series than the 

^ The summation in (2) runs over nonnegative integer j that are less than or equal 
to (j) + l){m — k) /p 
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number of numerical parameters in the Laurent series plus the number 

of hj k. 

(3) Choose a Laurent series expansion for some function and substitute the 
obtained Laurent series coefficients into Eq. (2). This substitution trans- 
forms (2) into a linear and overdetermined system in hj i. with coefficients 
depending on numerical parameters. 

(4) Eliminate coefficients hj ^, and get a nonlinear system in parameters. 

(5) Solve the obtained nonlinear system. 

R. Conte and M. Musette note that a computer algebra package is highly 
recommended for using of their method [9]. Steps of this algorithm can be 
implemented in computer algebra systems separately. 

For the given system it is easy to calculate the Laurent series solutions to 
any accuracy. These computations base on the Painleve test, which has been 
implemented in the most popular computer algebra systems [4,34,38,50]. Note 
that when one has computed a sufficient number of the Laurent series coeffi- 
cients he can forget about the system of differential equations and work only 
with coefficients of the obtained series. The first package of computer algebra 
procedures, which realize the third and the fourth steps of the algorithm, has 
been written in AMP [13] by R. Conte. One can also use our Maple and RE- 
DUCE packages of procedures, which are accessible in Internet [44] and are 
described in [45,46]. So, one passes the first four steps of algorithm without 
any difficulties. 

At the fifth (last) step one should solve an overdetermined system of nonlinear 
algebraic equations. The standard method for solving of such systems is the 
construction of a lexicographically ordered Grobner basis [11]. The Buchberger 
algorithm [5,6,11] guarantees a construction of a Grobner basis in a finite 
number of steps ^ . However, appropriate estimates of the number of steps 
required for constructing of this basis do not exist. The required computer 
memory depends, in the general case, exponentially on the number of the 
unknowns. Therefore, this number should be made as small as possible. 

The purpose of this paper is to show that we can essentially simplify the 
algebraic system of equations, which we have to solve on the last step of the 
Conte-Musette method, if we search the elliptic solutions only. 

In [32] R. Conte and M. Musette have used their method to find wave solutions 
of the complex cubic Ginzburg-Landau equation (CGLE3). The nonexistence 
of elliptic travelling and standing wave solutions of the CGLE3 has been proved 
in [23] and [47] respectively. In Section 4 we seek the eUiptic solutions of the 
quintic complex Ginzburg-Landau equation (CGLE5) using our modification 

^ An alternative method for constructing Grobner bases is the involutive algo- 
rithm [18,19]. 
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of the Conte-Musettc method. Note that both the CGLE3 and the CGLE5 
have only one-parameter Laurent-series solutions in the wave form and there 
exist only a finite number of such solutions. 

Our approach can be effectively used in investigation of any system (1) or a 
single differential equation, for which only a finite number of different Laurent- 
series solutions exist. Note that for wide class of such differential equations 
it has been proved that all their meromorphic solutions are elliptic (maybe 
degenerated) functions [15]. 



3 Properties of the elliptic functions 

Let us recall some definitions and theorems. The function q{z) of the complex 
variable z is a doubly-periodic function if there exist two numbers cui and 002 
with ojxjijji ^ such that for all 2; e C 

q{z) = g{z + cji) = g{z + cj2)- (3) 

By definition a double-periodic meromorphic function is called an elliptic func- 
tion [14,22]. These periods define the period parallelograms with vertices Zq, 
Zq + Nicui, Zq -\- N2L1J2 and Zq -\- NiUJi -\- N2UJ2, where A^i and N2 are arbitrary 
natural numbers and zq is an arbitrary complex number. The fundamental 
parallelogram of periods is called a parallelogram of period, which does not 
include other parallelogram of periods, that corresponds to A^i = A^2 = 1- 

The classical theorems for eUiptic functions [14,22] prove that 

• If an elliptic function has no poles then it is a constant. 

• The number of elliptic function poles within any finite period parallelogram 
is finite. 

• The sum of residues within any finite period parallelogram is equal to zero 
(the residue theorem). 

• If (){z) is an elliptic function then any rational function of g{z) and its 
derivatives is an elliptic function as well. 

• For each elliptic function g{z) there exist such m (m ^ 2) and such coeffi- 
cients hi J that g{z) is a solution of Eq. (2). 

Lemma 1 An elliptic function can not have two poles with the same Laurent 
series expansions in its fundamental parallelogram of periods. 

Proof. 

Let some elliptic function g{^) has two poles in points ^0 and ^1, which belong 
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to the fundamental parallelogram of periods. The corresponding Laurent series 
are the same and have the convergence radius R. Then the function v{^) = 
QiC ~ Co) ~ Qii ~ Ci) is an elliptic function as a difference between two elliptic 
functions with the same periods. At the same time for all ^ such that |^| < it! 
v{^) — 0, therefore, v{^) = and — ^o) = — Ci) Ci — Co is a period 
of g{C)- It contradicts to our assumption that both points and belong to 
the fundamental parallelogram of periods. 



4 Construction of elliptic solutions 

4-1 The quintic complex Ginzburg-Landau equation 

The one- dimensional quintic complex Ginzburg-Landau equation (CGLE5) 
is a generalization of the one-dimensional cubic complex Ginzburg-Landau 
equation [20] (CGLE3), which is one of the most-studied nonlinear equations 
(see [3] and references therein). Moreover, the CGLE5 is a generic equation 
which describes many physical phenomena, for example, the behaviour of trav- 
elling patterns in binary fluid convection [36] and the large-scale behavior of 
many nonequilibrium pattern-forming systems [10]. 

The CGLE5 is as foUows 

iAt + pAxx + q\A\'^A + r\A\^A - i-fA = 0, (4) 

where subscribes denote partial derivatives: At = Axx = ^r, P-,(l-,'<^ G C 
and 7 G M. 

One of the most important directions in the study of the CGLE5 is the con- 
sideration of its travelling wave reduction [2,9,12,31,35,36]: 

A{x,t) = ^M{i) e^i'PiO - ^i) ^ ^^^_ct^ ceM, ueR. (5) 

Substituting (5) in (4) and multiplying both sides of this equation on AM'^/A 
we obtain 

2pM"M - pM'^ + AipnpMM' + 2(2t^ - ic - 2i7 + 

(6) 

+ 2cV' - 2pV^ + 2ipV')M2 + AqM^ + ArM"^ = 0, 

where ijj = ip' = M' = Equation (6) is a system of two equations: 
both real and imaginary parts of its left-hand side have to be equal to zero. 
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Dividing (6) by p and separating real terms from imaginary ones, we obtain 
the following system 

2MM" - M'^ - AM^ip^ - 2cMM' + Ag^M^ + Ad^M^ + 4m^M^ = 0, 

(7) 

MV-' + -0 (M' - cM) - QrM + diM'^ + UiM^ = 0, 



where new real variables are as follows 

r q \ 

Ur-\-mi — -^ dr + idi — -, Sr — iSi — -, (8) 

p p p 

9r + m = (7 + i^) {Sr - iSi) + ^C^SiSr + |C^S^, (9) 



and 



ijj = - g = csi. (10) 



System (7) includes seven numerical parameters: Qr, gi, dr, di, Ur, Ui and c. 
Note that to obtain (7) from (6) we assume that the functions M(^) and 4>{^) 
are real. 

The standard way to construct exact solutions for system (7) is to transform 
it into the equivalent third order differential equation for M. We rewrite the 
first equation of system (7) as 

V^' = (11) 



where 



G = -MM" - - -MM' + giM^ + drM^ + UrM^. (12) 



Prom (11) it follows that 
~,~ G'M-2GM' 

V"0 = 5 , (13) 

2M3 ' ^ ^ 



Multiplying the second equation of (7) on ip and substituting (11) and (13) in 
it, we express ip in terms of M and its derivatives: 

J _ G'-2cG 

^ ~ 2M2 {g, - diM - UiM^) ' ^ 
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and obtain the third order equation for M: 

{C - 2cGf + AGM^{gr - d,M - u.M^f = 0. 



(15) 



4-2 The Laurent series solutions 



Below we consider the case 

-^K, (16) 
r 



which corresponds to the condition Ui ^ 0. In this case Eq. (15) is not inte- 
grable [31] and its general solution (which should depend on three arbitrary in- 
tegration constants) is not known. Using the Painleve analysis [31] it has been 
shown that single-valued solutions of (7) can depend on only one arbitrary 
parameter. System (7) is autonomous, so this parameter is ^q: if M = /(^) is 
a solution, then M = f{C, — C,o), where C,o ^ has to be a solution. Special 
solutions in terms of elementary functions have been found in [2,9,31,35]. All 
known exact solutions of (7) are elementary (rational or hyperbolic) functions. 
The full hst of these solutions is presented in [9] . The purpose of this section 
is to find at least one elliptic solution of (7). 

System (7) is invariant under the transformation: 

■0 ^ - -0, Qr^ - Qr, - di, Ui^ - Ui, (17) 



therefore we can assume that ui > Q without loss of generahty. Moreover, 
using scale transformations: 

M AM, dr ^ ^, di^ ^, Ur ^7^, ^, (18) 



we can always put ui — l. 

Let us construct the Laurent series solutions to system (7). We assume that 
in a sufficiently small neighborhood of the singularity point the functions 
■0 and M tend to infinity as some powers of ^ — ^o- 

^ = A{i-ioT and M = i?(e-eo)^ (19) 



where a and j3 are negative integer numbers and, of course, A and B ^ 0. 
Substituting (19) into (7) we obtain that two or more terms in the equations 
of system (7) balance if and only if a = —1 and (3 — —1. In other words in 
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this case these terms have equal powers and the other terms can be ignored 
as t — > tg. We obtain values of A and B from the following algebraic system: 

52 (3 - AA^ + 4urB^) = 0, 

(20) 

2A-B^^ 0. 



System (20) has four nonzero solutions: 



Ai^Ur + ^^4x^2 + 3, ^1 = V + + (21) 



A2^Ur + -^jAul + 3, = - \l2ur + Au^ + 3, (22) 



A^ = Ur-- ^|4.ul + 3, 53 = \l 2ur - ^jAul + 3 (23) 



and 



A^ = Ur- --^^ul + 3, 54 = - ^j2ur - ^Aul + 3. (24) 



Therefore, system (7) has four types of the Laurent series solutions. Denote 
them as follows: 

A B 
V^fe = y + «M + «M^ + ---' ^fc = y + ^M + M + ---' (25) 



where k = l.A. 

Let M(^) is a nontrivial elliptic function. Note that if ■?/' is a constant, then from 
the second equation of system (7) it follows that M can not be a nontrivial 
elliptic function. Therefore, using (14), we conclude that has to be a 
nontrivial elliptic function as well. 

Let us consider the fundamental parallelogram of periods for the function 
M(^) and define a number of its poles in this domain. Let M has a pole of 
type Ml, hence, according to the residue theorem, it should has a pole of type 
M2 (it can not have a pole of type M4 because Ur is a real parameter) . So ip 
has poles with the Laurent series ■0i and ■02- As an elliptic function it should 
have a pole of type ip^ or ip^ as well. It means that the function M{^) should 
have a pole of type M3 and, hence, a pole of type M4. So M{^) should have at 
least four different poles in its the fundamental parallelogram of periods. Using 
Lemma 1, we obtain that the function M(^) can not have the same poles in 
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the fundamental parallelogram of periods. Therefore, M(^) has exactly four 
poles in its fundamental parallelogram of periods. In this case by means of the 
residue theorem for ip we obtain 

Ur = 0. (26) 



We obtain that the CGLE5 with 7^ has no elliptic solution in the wave 
form. In the case Ur — Q possible elliptic solutions should have four simple 
poles in the fundamental parallelogram of periods, and, therefore, has the 
following form [22]: 

M{i-i^)^c + j2Bkai-ik). (27) 

k=l 



where the function C(^) is an integral of the Weierstrass elliptic function mul- 
tiplied by —1: 

C'(0 = -P(0, (28) 



C and are constants to be defined. We also should define periods of the 
Weierstrass eUiptic function. 

To obtain restrictions on other parameters, we use the Hone's method [23] 
and apply the residue theorem to the functions and so on. The residue 

theorem for the function t/;^ gives the equation: 

fe=i 



The values of follows {Ur — 0): 



«1,0 



^4,0 



^ (6c - </27di - 154^dr) , (30) 



02 = ^ (^6c + </27di + 15^0?,) , (31) 



ag^o = - + i {</^di - 15</3dr)) , (32) 



-^{6c-i(</27di-15</3dr)) . (33) 
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Substituting and 0^,0 in (29), we obtain 

E^fc«M = 25 = 0, (34) 

k=l ^ 



therefore c = 0. 

For the function ijj^ the residue theorem gives: 

4 

^^fe(Aafe,i + 4_o) =0. (35) 

k=l 

The values a^. ^ are as follows (we put c — 0): 

^11 = (3o?j^ — 63dr^ — ddVSdidr + 128V3gi) , ^2 1 = ^1 (36) 

' 384 ^ / ' ' 

03,1 = (Sdi^ - 63dr^ + QdVSdidr - 128V3gi^ , a^ ^ = a^ ^. (37) 



Equation (35) is equivalent to 

di^ + 27dr^ = ^ di = ±iV27dr. (38) 

The parameters dr and rfj should be real, therefore, dr — and di — 0. So, 
consideration of ijj'^ and gives three restrictions: 

c = 0, 4 = and di = 0. (39) 

The residue theorem for ijj^ gives the restriction 

QiQr = 0. (40) 

Considering ip^ and ijj^, we do not obtain new restrictions on coefficients. 
Taking into account (26) and (39) we obtain system (7) in the following form: 

2MM" - M'^ - AM^"^ + AgiM^ = 0, 

(41) 

ip'M + ipM' - QrM + = 0. 

To find elliptic solutions to system (41) we use the Conte-Musette method. 
The function ip{^) can have not four but two different Laurent series expan- 
sions, whereas the functions M(^) should have four different Laurent series 
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expansions, so it is easier to find than Af(^). Equation (2) with m — 1 
has no eUiptic solution. Let satisfies Eq. (2) with m — 2: 

V^'^/ 



Substituting in (42) the Laurent series of ip, which begins from Ai (more 
exactly we use the first ten coefficients), we obtain the following solution hkj 
for an arbitrary value of the parameter Qr ^ and gi — 0: 

4 _ ^2 ^ . . . . . 

^4,0 — ~ i:, ho,o — — /i3,o = ^2,0 = ^1,0 = ^0,1 = ^1,1 = ^2,1 = 0, (43) 



a few solutions with Qi — and Qr — and no solution for gi 0. 

The straightforward substitution of the functions 



or 



2t ' 



M 



±i 



(45) 



in (41) with Qr = and Qi = proves that they are exact solutions. The 
coefficients of the Laurent-series solutions do not include arbitrary parameters, 
so the obtained solutions are unique single- valued solutions and system (41) 
has no elliptic solution for these values of parameters. 

In the case of solutions (43) the function satisfies the equation 



The polynomial in the right hand side of (46) has four different roots, therefore 
ip is a non-degenerate elliptic function [22]. 

Surely we do not rigorously prove the existence of elliptic solutions to the 
CGLE5. More precisely, we calculate fifty coefficients of the Laurent series of 
the function ip, which satisfies (41) with Qi = and find that they coincide 
with corresponding coefficients of the Laurent series of the exact solution to 
Eq. (46). 
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For rigorous proof we should find the function M(^) and check that this func- 
tion is a sohition of (15). The function Af(^) in a parallelogram of periods has 
four different Laurent series expansions, so we should choose the parameter m 
such that solutions of Eq. (2) have four poles in its fundamental parallelogram 
of periods. Minimal possible value of m is equal to 4. The general form of (2) 
for m = 4 and p = 1 is the following: 



Substituting the Laurent series Mk from (25), we transform the left hand side 
of (47) into the Laurent series, which has to be equal to zero. Therefore, we 
obtain the algebraic system in hij and gr- The first algebraic equation, which 
corresponds to 1/^^ is 



where is defined by (21)-(24). If we would attempt to find all elliptic 
and degenerate elliptic solutions then we should use only one of Bj. and can 
express, for example. /?.8 q via /ig i, /i4 2 and /i2.3. We seek only elliptic solutions, 
so we know that all B^. have to satisfy (48) and can consider Eq. (48) as the 
following system: 



M'^ + (/l2,3M2 + /li,3M + ho,3) M'^ + 

+ {h4,2M^ + /l3,2M3 + h2,2M^ + + /lo,2) + 

+ (/l6,l^^ + Kl^'" + K^iM^ + /l3,iM3 + /l2,lM2 + hi^iM + 

+ /io,l) M' + /i8,oM8 + h-j^oM'' + /i6,oM6 + /i5_oM^ + 

+ /i4,oM^ + /i3,oM3 + /i2,oM2 + Hi^qM + /iq.o = 0. 



(47) 



Bl {hsflBl - he,,Bl + h,,2Bl - h2,3Bk + l) = 0, 



(48) 



< 




he,iBf + K,2Bl 

/l6,l-B| + hi2B2 

h&,\Bl + /i4,2-B| 



h2,3Bi + 1=0, 
h2,3B2 + 1=0, 
h2,3B3 + 1 = 0, 
/l2,354 + 1=0. 



(49) 



Using the exphcit values of B^ from (21)-(24), we obtain that 



^8,0 — — o' ^4,2 



= 0, he,i = 0, /i2,3 = 0. 



(50) 
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Taking into account (50), from other equations of the algebraic system we 
obtain 

J, _ 4 , _ 16 2 ^ _ 64 3 

1^6,0 — -^9r, ii4,0 — '^2,0 — l^-'-J 



all other hij are equal to zero. So, the equation for M has the form 

M'^ = —M^ (SM^ - Ag^y . (52) 



Equation (15) at Ui = 1, Ur = 0, c = 0, dr = 0, di — and Qi — has the 
form: 

1 (^M"'f - (2MM" - M'^) (m^ -grf ^ 0. (53) 



We multiply Eq. (53) on M'^ and use Eq. (52) to express all derivatives of 
M(^) in terms of the function A'/(^). The straightforward calculation shows 
that any solution of (52) satisfies (53). So, we obtain elliptic wave solutions of 
the CGLE5. If Sj 7^ these solutions are the standing wave solutions, in the 
opposite case (sj = 0) the solutions can have an arbitrary speed c. 

Note that we obtain (15) from (6) using the condition that M(^) is a real 
function. For gr < and any initial value of M we obtain real solutions. In 
the case gr > there exists the minimal possible initial value of M for which 
real solutions exist and only particular solutions of (52) are suitable elliptic 
solutions to the CGLE5. The function M(^) has the form (27), the values of 
constants can be determined from (52). 

Summing up we can conclude that our modification of the Conte-Musctte 
method allows us to get two results: we obtain new elliptic wave solutions of 
the CGLE5, and we prove that these solutions are unique elliptic solutions for 
the CGLE5 with gr 7^ 0. 

From (40) it follows that elliptic solutions can exist if gr or gi is equal to zero. 
For all nonzero values of gr and zero gi we have found elliptic solutions. In 
the case gr — 9% — ^ there is no elliptic solution. In the case of zero gr and 
nonzero g^ we substitute the obtained Laurent series solutions M^. into Eq. (2) 
with 171 = 1, . . . , 4 and do not obtain neither elliptic functions nor degenerate 
elliptic solutions. We hope that the more detail analysis of this case allows us 
to find all eUiptic solutions for the GGLE5. 
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5 Construct of elliptic solutions for nonintegrable systems 

The approach, which we have considered in the previous section, can be apph- 
cable to many nonintegrable systems. The best applicability is nonintegrable 
systems of autonomous nonlinear ordinary differential equations with so-called 
finiteness property [15]: there is only a finite number of formal Laiuent series 
that satisfy the system. For such systems we can propose the following way 
to the search for eUiptic solutions: 

(1) Calculate a few first terms of all solutions of system (1) in the form of 

the Laurent series. 

(2) Choose the function t/^, which should be elliptic. Check should other 
functions be elliptic or not. 

(3) Using the residue theorem define values of numeric parameters at which 
the solution yk can be an elliptic function. 

(4) Define a minimal number of poles for candidates to elliptic solutions. 
Using this number, choose a positive integer m, define the form of Eq. (2) 
and calculate the number of unknown coefficients hj^]^. 

(5) Calculate the sufficient number of coefficients for all Laurent series of 
Hk and substitute the obtained coefficients into Eq. (2). This substitu- 
tion transforms (2) into a linear and over deter mined system in with 
coefficients depending on numerical parameters. 

(6) Eliminate coefficients ^ and get a nonlinear system in parameters. 

(7) Solve the obtained nonlinear system. 

We restrict ourselves to the search of elliptic solutions only. This makes possi- 
ble to consider not one Laurent series, but as much as possible. This is a key 
idea of our modification of the Conte-Musette method. 

Using the Conte-Musette method in the original form, one in principal would 
be able to find the elliptic solutions of the CGLE5, but calculations would 
be cumbersome, because he should reobtain all known solutions in terms of 
elementary functions to obtain new elliptic solutions. One has to construct and 
to solve an algebraic system, which includes at least 32 equations, which are 
linear in 24 unknowns hi^ and nonlinear in 7 parameters. In our approach, we 
seek only elliptic solutions and, using the Hone's method, reduce the number 
of arbitrary numerical parameters to 1. Consideration of a few Laurent series, 
instead of one also simplifies the obtained system of algebraic equations. We 
conclude that the use of the Conte-Musette method in combination with the 
Hone's method is very effective. 

In contrast to [23,32] we consider a system of differential equations instead 
of the equivalent single differential equation and demonstrate that the system 
is more convenient for analysis than a single equation for M. If we analyse 
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only the Laurent series expansions of the function M, then we should also 
consider the case when M has not four, but only two different Laurent series 
expansions, beginning from Bi and B2 or and B^. The consideration of the 
Laurent series expansions of two functions: and M, allows us to reject this 
possibility. Note that in our opinion if one seek degenerate elliptic solutions 
using the Laurent series then a system of equations is more convenient for 
analysis than the equivalent single equation as well. 

Moreover, in contrast to traditional methods when we use the Conte-Musette 
method we can choose a function, which analytic form should be found. Instead 
of the functions ip and M we can consider some combination, for example, a 
polynomial, of this functions and seek this combination in the analytic form. 
Note that we have no need of any differential equation for this combination. So, 
we can conclude that the use of the Laurent series solutions gives new insight 
on problem of construction of exact solutions for nonintegrable systems. 



6 Compcirison with standard methods 

Let us compare our approach with other methods for construction of special 
solutions of nonintegrable systems (see [9] as a review of these methods) . 

Without the use of the Laurent series solutions it would be very difficult to 
find the elliptic solutions of the CGLE5, because the form of Eq. (47) is very 
complex. For example, the Fan's technology [16,17] allows to find solutions of 
nonintegrable systems, which are polynomials of the function q, that satisfies 
the following equation: 

TV 

/ = ^c,^, (54) 
3=0 

where N and Cj are constants to be determined. The function ip is a. solu- 
tion of a similar equation with N — 4, but, using the Fan's approach, it is 
impossible to find without knowledge of M(^), which satisfies the more 
complex equation (47). The Kudryashov's method [26] as well as the meth- 
ods [9,16,17,21,25,27,32,33,37,40,41,43], proposed to search both elliptic and 
elementary solutions, attempt to find solutions for the initial differential sys- 
tem at all values of the numeric parameters such that single- valued solutions 
can exist. The use of the residue theorem and the Hone's method allows us 
to fix some of these parameters and to simplify calculations without loss of 
elliptic solutions. Note that we not only find elliptic solutions for the CGLE5, 
but also we prove that there are no other elliptic solutions for ^ Using 
the standard methods one can say nothing about uniqueness of the obtained 
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special elliptic solutions. 



To find solutions for the CGLE5 the authors of papers [2,31,35] put some 
restrictions on dependence between ip and M. The use of Laurent solutions 
allows us to search ijj without any restrictions and without eliminating M from 
system (7). 



7 Conclusions 



In this paper we propose a new approach for the search of elliptic solutions to 
systems of differential equations. The proposed algorithm is a modification of 
the Conte-Musette method [32] . We restrict ourselves to the search of eUiptic 
solutions only. A key idea of this restriction is to simphfy calculations by 
means of the use of a few Laurent series solutions instead of one and the use 
of the residue theorem. 

The application of our approach to the quintic complex one-dimensional Ginz- 
burg-Landau equation (CGLE5) allows to find elliptic solutions in the wave 
form. Let us point out that the obtained solutions are the first elliptic solu- 
tions for the CGLE5. Using the investigation of the CGLE5 as an example, we 
demonstrate that to find elliptic solutions the analysis of a system of differ- 
ential equations is more preferable than the analysis of the equivalent single 
differential equation. 

We also find restrictions on coefficients, which are necessary conditions for 
the existence of eUiptic solutions for the CGLE5. To do this we develop the 
Hone's method [23] . We show that this method is useful not only to prove the 
nonexistence of elliptic solutions, but also to find new elliptic solutions. Note 
that the Hone's method and, therefore, our approach, are so effective in the 
case of the CGLE5, because coefficients of the Laurent series solutions depend 
only on parameters of equations, i.e. they do not include additional arbitrary 
parameters (have no resonances). It is an important problem to generalize the 
Hone's method on the Laurent series solutions with resonances. 



Another way for future investigations is the generalization of the Conte- 
Musette method on the case of multivalued solutions. Some results in this 
direction have been obtained in [43,48]. 
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